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Abstract 

L». ' In this work we show that 3d Feynman ampHtudes of standard QFT in flat and homogeneous space can 

^— ^ , be naturally expressed as expectation values of a specific topological spin foam model. The main interest of 

^SJ ' the paper is to set up a framework which gives a background independent perspective on usual field theories 

» , , and can also be applied in higher dimensions. We also show that this Feynman graph spin foam model, 

C^ ' which encodes the geometry of flat space-time, can be purely expressed in terms of algebraic data associated 

with the Poincare group. This spin foam model turns out to be the spin foam quantization of a BF theory 

based on the Poincare group, and as such is related to a quantization of 3d gravity in the limit Gjv —^ 0. We 

investigate the 4d case in a companion paper llj where the strategy proposed here leads to similar results. 
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I. INTRODUCTION 

There has been a recent understanding and new results concerning the consistent couphng of 
matter fields to 3d quantum gravity amplitudes, in the context of the spin foam approach to 3d 
quantum gravity [^, y, \4, la, la, [3] . These results have led to the construction of an effective field 
theory, describing the coupling of quantum fields to quantum gravity [3, ISJ , which arises after the 
exact integration of quantum gravitational degrees of freedom. This effective field theory lives on 



a non-commutative space-time whose geometry is studied in detail in [10(]. When the no gravity 
limit Gn ^ is taken, the geometry becomes commutative and usual field theory is recovered. 

Prom the spin foam side, the way usual Feynman diagrams are recovered is highly non trivial, 
since quantum gravity amplitudes are expressed as purely algebraic objects in terms of spin foam 
models, a language quite remote from usual field theory language. The fact that Feynman diagrams 
arise from this picture effectively, amounts to show that usual field theory can be expressed in a 
background independent manner, where fiat space-time emerges dynamically from the choice of 
spin foam amplitudes. Namely, this amounts to show that Feynman amplitudes can be expressed 
as the expectation value of certain observables in a topological spin foam model. 

Note that this last statement is not a statement regarding quantum gravity but one regarding 
well known Feynman diagrams in 3d. It is then natural to wonder if one can reach a deeper under- 
standing of this property; whether it is tied up to the topological nature of gravity in dimension 3 
or if it can be extended to field theory in higher dimensions. It also begs the question, whether or 
not Feynman diagrams contain some seeds of information about quantum gravity dynamics. 

A natural strategy to address successfully these questions is to start from a careful study of 
the structure of Feynman graph amplitudes, without any assumption about the nature of quantum 
gravity dynamics. The goal is to see why and how such familiar objects can be interpreted in terms 
of spin foam models, which have appeared to be a suitable language to address dynamical question 



in a background independent approach to quantum gravity [11|. 

Usual Feynman diagrams are natural observables^ from which physical predictions regarding 
matter interactions in space-time can be extracted. As such, they allow one in principle to probe 
the geometry of space-time at arbitrarily small scales. Moreover, they can be naturally promoted 
to physical observables in any theory of quantum gravity. Therefore, the preferred spin foam 
model - if any - which is related to Feynman diagrams should give us some important seeds of 
information about the structure of quantum gravity amplitudes. In other words, if one assumes 
that it is possible to consistently describe quantum gravity amplitudes coupled to matter in terms of 
a spin foam model, a necessary requirement for the spin foam amplitudes, in order to be physically 
relevant, is to reduce to the ones appearing in the study of Feynman diagrams when Gn — > 0. This 
gives strong restrictions on the admissible, physically viable spin foam models. 

In this paper we focus our study to the case of 3-dimensional Feynman diagrams, but as should 
be clear from our work and as shown in a companion paper pj , most of the results obtained here 
can be extended to the more interesting and challenging case of 4-dimensional amplitudes. What 
makes the analysis work is not that gravity is a topological theory but the fact that, when coupled 
to matter, quantum gravity in the limit Gm — > behaves as a topological field theory. This is a 
trivial fact in dimension 3 and has been argued to be true in dimension 4 when one takes this limit 
while preserving diffeomorphism invariance [12l |. 

The main idea behind our analysis is the fact that in 3d, quantum gravity amplitudes provide an 



^ Of course, any observable quantity is actually given by a sum of Feynman diagrams, for instance those of fixed 
valency and given degree, as generated by a field theory at given order of perturbation, and subject to renor- 
malization procedure. In order to keep the exposition simple we will refer to individual Feynman diagrams, in a 
regularized form in order to avoid divergences. 



integration measure for the integration of Feynman graph amphtudes. This idea is contained in [3] 
and spelled out in [l3|. When gravity is turned off, this measure should be related to the Lebesgue 
measure expressed in terms of relatives lengths, which are Poincare scalars. When gravity is turned 
on, this measure is deformed and reveals the quantum nature of the corresponding space-time. 

The organization of this paper is as follows. In section |TI] we first express explicitly the change 
of variables from the Lebesgue measure that arises in Feynman graph integration, to the measure 
expressed in terms of relative distances. In section IIIII we present our main result, namely that 
Feynman integrals can be expressed as the expectation value of an observable for a topological 
spin foam model. This model is constructed out by removing all explicit information about flat 
space into a uniquely defined choice of amplitudes. We carefully study this model, especially its 
symmetries and its gauge fixing properties, and show that it defines a topological state sum model. 
This topological model, which naturally emerges from Feynman graphs, has already been studied 



in the mathematical literature from a different perspective by Korepanov IJ] . In section IIVI we 
show that the tetrahedral weights entering the definition of the topological state sum model can be 
given an algebraic interpretation in terms of the 6j-symbol of the Poincare group; this 6j-symbol is 
computed explicitly and for the first time. We also relate this model to the Ponzano-Regge model 
via a doubling and limiting procedure. Finally we show how to extend these results to the case of 
a non vanishing cosmological constant. 

II. FEYNMAN DIAGRAMS IN TERMS OF INVARIANT MEASURE 

We consider QFT in flat Euclidean D-dimensional space-time and focus on closed Feynman 
diagrams for a scalar field. The Feynman amplitude for a graph F is given by an integral over the 
position in M of N vertices. The (Euclidean) Poincare group ISOd acts diagonally on these A'^ 
points 

(xi , • • • ,xn) ^ {A- xi + a,- ■ ■ ,A- XN + a) 

where (A, a) G SOd x ^^ is a general element of the Euclidean Poincare group. The integrand is 
a function of the distances between the vertices, namely an invariant under the action of ISOd 

It= I d^xi---dPxNOT{\xi-Xj\), with Or = n G^{xi-Xj) (1) 

(u)Gr 

Gp is the Feynman propagator and the product is over all edges of F. Following the ideas of 
[a, ll3] , we want to express this integral in terms of the invariant measure, involving the quantities 
lij = \xi — Xj\. On one hand, ND variables are integrated over in ([1]) while the dimension of the 
Poincare group is ^ D{D + 1); consequently, the invariant measure depends on ND — ^ D{D + 1) 
variables: the number of independent variables required to place the vertices with respect to each 
other^. On the other hand, the number of length variables lij is ^AfiN — l). Note that the quantity: 



N{N - 1) 



iVP-°'° + " 



(iV - 0)(JV - P - 1) pj 



which computes the number of redundant edge lengths, vanishes only for N = D and N = D + 1. 



^ Precisely this quantity is the number of continuous variables required to place the vertices with respect to each 
other modulo discrete transformations. 



First, for N = D, we split the measure into the Haar measure for ISOo and a product of Ujdlij 

d^xi ■ ■ ■ d^XD = d^a dA JJ kjdkj, (3) 

i<j 

where a is the barycenter of the D points and dA is the Haar measure on SOd- The volume of 
SOd with respect to this measure is given by 

Ad 



ntiTd) ' 

so A2 = 2tt, A3 = 2^7r^ and A4 = 2 vr . Given an additional point xd+i, one can build a /^-simplex 
in M^ with vertices Xj, and express the Lebesgue measure in terms of the edge lengths as 

,D _ V^ Ui=l liD+ldljD+l , . 

°"",k — ^) — " 

V{lij) is D! times the volume of the simplex whose edge lengths are lij. In the rest of the paper, 
we will simply call 'volume' the quantity V. e labels the orientation'^ of the D-simplex xi- ■ ■ xd+i- 
The proof of (j3]) is easily obtained by computing the jacobian of the linear function which maps 
an orthonormal basis Cj to the basis defined by the vectors {xi}i=i...£)- Let y G M and li = \y — Xi\; 
we have 

D 



= n^^* |det(ej • {y - Xi))\ = d^y^Jdet{{y - Xj) ■ {y - Xi)) (6) 

This last determinant is Dl times the volume of the simplex xi, • • • , XD,y- The previous identity is 
a local identity showing that the map y — > ij is locally invertible away from the hyperplane spanned 
by xi, • • • , XD- This map is globally invertible on each half space separated by this hyperplane. The 
sum over orientations comes from the fact that the set of edge lengths li determines the position 
of y up to a reflection with respect to the hyperplane spanned by xi, • • • ,xd- Such a reflection 
is in 0{D), but not in SO{D), and changes the orientation of the D-simplex. The equality of 
measures (JH) is then understood as an equality for integrating out a function of x/5+1 in the LHS 
and a function of hn+i and orientation e in the RHS. The proof of ([3]) is obtained by a direct 
computation and by recurrence. 

So putting (j3|4p together we get, for N = D + 1 points in M , that the invariant measure is 
given by: 

d xi---d XD+i =d adA 2^ ' . (7) 



A similar identity already appeared in 15 1. 



Now for N = D + 2, the quantity ([2]) is equal to 1. Accordingly, the invariant measure will be 
formulated using all the distances, except for one. One can indeed express it in terms of the edge 
lengths of the polyhedron formed by two D-simplices glued together in M along a {D — l)-simplex 



The notion of orientation is precisely defined in appendix 1X1 



(as shown in Fig [2] for D = 3). The missing edge length is the one joining the only two vertices, 
denoted by a,b, which do not belong to the common (D — l)-face. The measure is obtained by 
using twice (JH) 



d^xi- 



• d^XZ)+2 



d^adA 









(8) 



where the product is over all edges {ij) different from (a, b), Va denotes the volume of the Z'-simplex 
not containing a, and the sum is over orientations of the D-simplices a, b. 

Note that the edge lengths kj, (ij) ^ (ab) only determine the geometry of the polyhedron 
modulo reflections of the D-simplices along the common (D — l)-face. In order to fully determine 
the geometry, one needs to specify the orientations of the two D-simplices denoted a and b. Once 
these orientations ea, e^ are specified, the missing edge length can be uniquely reconstructed: lab = 
^afe(^jj,e)- For fixed kj, the missing distance lab can take two values denoted by l^j^, with /^^ < /^^, 
depending on whether the two points (o, b) are separated or not by the hyperplane spanned by 
the common face. A simultaneous change of the two orientations is equivalent to a symmetry with 
respect to this hyperplane, under which all distances remain invariant. Therefore, the length lab 
depends only on the product eaCb, in such a way, according to our conventions, that l^^^ = I'al'^''- The 
Fig [1] illustrates the D = 2 case. More details are given in appendix. 





FIG. 1: Two triangles acd and bed in the Euclidean plane share an edge (cd). 



On the left (resp. right) 



For a general value of A^ = D + k, we use repeatedly {k times) the equation ^ in order to 
get the invariant measure. This invariant measure is based on a triangulation which is constructed 
recursively: each time a vertex is added, we choose a {D — l)-face on the existing triangulation, 
and build a new i^-simplex from this data. At each step we use ([4]) to construct the invariant 
measure whose form depends on the triangulation A^. The triangulations obtained this way are 
of special type: they are triangulations with a boundary being topologically a {D — l)-sphere and 
such that every {D — 2)-face lies on the boundary. 

More formally, let B^ be a Z)-dimensional ball B^ with n = D + k points marked on its 
boundary. We now consider a triangulation A^ of Bo with D + k vertices, such that every vertex 
and every {D — 2)-face are on the boundary. In other words A^ is a triangulation of B^ which 
possesses D+k boundary vertices and which does not possess any internal (D— 2)-face. If one draws 
the dual 2-skeleton of A^, in which vertices are dual to D-simplices, edges to {D — l)-simplices 
and faces to {D — 2)-simplices, one sees that our condition on A^ implies that this dual graph is a 
D-tree (a D + 1 valent graph containing no closed loop) with open ends. It can be checked that 
the number of edges of such a triangulation matches the number of Poincare invariant quantities. 



Given such a triangulation the invariant measure reads 

d^xi • • • d'^xD+k = d^adA Y. n '"^^^^ n ^ ^^^ 

ee{±l}'=eGAfe creAft " 

where the first product is over all edges e of the triangulation A^, the second product is over all 
its D-simplices a and the summation is over all orientations of its D-simplices. 

By construction the data {/ei £0-} defines uniquely a flat geometry on A^ and allow us to recon- 
struct the 'missing' edge lengths, namely the length of the edges that do not belong to A^. The set 
of distances {Uj} in the integrand of ([1]) splits up into the lengths {le-,e € A^,} and the 'missing' 
distances denoted by {l\i}- These missing distances are e-dependent functions of the le- 

The Feynman amplitude ([1]) is flnally given by: 

IV= I Wledle Y. U ^(^rileJl^ile)) (10) 

where we have dropped out the overall factor which corresponds to the gauge volume. 

Let us conclude this part: we have started with a Feynman amplitude on flat space-time; in 
the usual formulation, information about flat geometry is encoded in the Lebesgue measure. By 
working with the invariant measure, we have moved this information in the weight assigned to 
the triangulation, and in the relations lli{le) between 'missing' lengths, labels and orientations. 
Our goal now is to show that this information about flat geometry can be encoded entirely into a 
dynamical topological state sum (spin foam) model. 

III. FEYNMAN DIAGRAMS AS SPIN FOAM MODELS 
A. The main result 

One of the main results of this paper is the fact that the Feynman amplitude (jlOp can be written 
as the expectation value of an observable 



eer 
for the spin foam model: 



(11) 



^ ^ •' eeA {Se.er} \ ^ ^ / 

A is an arbitrary triangulation of a closed 3D-manifold and F is embedded into the one-skeleton 
of A. |e| is the number of edges of A. The edges carry a double label {le,Se), where le is a real 
positive number and Se is an integer. These labels are summed over a domain where triangular 
inequalities are satisfied for the /g's- e,- = ±1 denotes the orientation of the tetrahedron r. The 
product in parenthesis is over all tetrahedra. The action for each tetrahedron is similar to the 
Regge action ^] 

Sr{Se,le) = Y'^-^eile), (13) 

eGr 

with 61 (le) the interior dihedral angle of the edge e in r. The observable is simply the product of 
propagators in which the distances are replaced by the labels le living on the graph F. The equality 



pT]) is obtained if one restricts to trivial topologies (A is a triangulation of the 3-spliere S^). Note 
that all edges are summed over, and therefore no reference is made to flat geometry. 

Furthermore, we will show that this model is topological and that the weight associated to it 
can be interpreted algebraically as the 6j-symbol of the Poincare group. The rest of the paper 
is devoted to the demonstration of these statements. The notations we use are summarized and 
explained in details in appendix lAl 

B. Emergence of the model 

1. A key identity 

The keystone of the result announced above is a geometric identity associated with the possibility 



to embed a 4-simplex in M^. A similar identity appears in lj,ll7| as a basis to construct 3D manifold 
invariants. Let G = V^ be the square of the volume of a 4-simplex (0, . . . , 4) with edge lengths 
lij . We denote by Tj the tetrahedron obtained by dropping the point j of the 4-simplex, and by Vj 
its volume. We will consider both the complex of the three tetrahedra ti,T2 and t^ surrounding 
the edge (04), and the complex of the two tetrahedra ro,r4 sharing the face [123] (see Fig[2j). The 
orientation of tj is denoted by Cj . 

We work with fixed lengths lij,{ij) 7^ (04), the length /o4 being free to fluctuate. Then, the 
following identity of measures holds: 



4/o4<5(G)-2^ — 2^ ^04 ,,,,,, • (14) 



V0V4 ^^^ ""^ V1V2V3 



The delta function in the third term is the 27r-periodic delta function. The quantity loq^^ considered 
as a function of /o4 and orientations, is the deficit angle carried by the edge (04), defined to be: 



Wi 



E^^-^04 



04 — /^ tj 1/04 

i=i 

where 9q^ G [0)'^] is the interior dihedral angle of the edge (04) in tj. It represents the curvature 
localized on the edge, in the sense of Regge calculus. It vanishes modulo'^ 2tt if and only if the 
complex of three tetrahedra (Ti,T2,r3) can be mapped in M^ with the orientations ej. Thus, the 
delta function of the deficit angle acts as a projector on the space of flat geometries. 

On the other hand, the complex of the two tetrahedra ro,r4 can always be mapped in R^, 
in a unique (modulo translations and rotations) way depending on the orientations. The 'miss- 
ing' distance between the points 0,4 is then a well defined function I'q^'^ of lij and the product 
orientations. 



We now give a proof of (I14p . We will need three intermediate results: the first one states that 
the equation 



G{k 



04 j 



■* If we don't work modulo 2-r the condition for the existence of an embedding of the 4-siniplex (0, . . . , 4) in flat 
space is ciJo4 = 27rn(e) where n(e) G {0, +1, —1} depends on the orientations in this way; 

n[e) = h '' '' = '^ ^^'^^^ (15) 

^ ^ 1 if not ^ ' 




FIG. 2: Pentagonal representation of a complex of two tetrahedra 0,4 sharing a face [123], and a complex 
of three tetrahedra 1,2,3 sharing an edge (04). 

where G = V^ is considered to be a function of Iq^, all the other lengths being fixed, admits exactly 
two solutions /q^. It was implicit in our discussion above, where we have mentioned, moreover, that 
these solutions are labelled precisely by the product eoe4. This result is geometrically clear since 
for each value of /o4 one can embed the simplex in M^ in such a way to map, say, tq into a fixed 
tetrahedron. The volume of the simplex vanishes only when the vertex belongs to the space E 
spanned by the tetrahedron tq. As /o4 varies within its domain, the vertex moves along a circle 
around the face [123], intersecting E at exactly two points. 

We can also give an algebraic justification for having only two solutions. Indeed, G is the Cayley- 
Menger determinant of the 4-simplex. The Cayley-Menger matrix is a 6 x 6 symmetric matrix such 
that its diagonal elements are and the off diagonal elements are given by Lab = "^afe/^' except 
for the first column and row where they are 1. Its determinant is therefore a quadratic function 
in Lo4. Moreover, the coefficients of this polynomial can be computed explicitly (see ()A13P and 
(lAT5l) ): 

-2VoV4Cos0O4, 7T7^ = -2Vo'4, (16) 



9Lo4 ' dLl^ 

where ^04 is the interior dihedral angle between the two tetrahedra and 4, and Vo4 is the area of 
the face [123]. 

The second result states that the equation 

'^04(^04) = 0, 

where lo is considered as a function of /o4 (the other lengths being fixed) and orientations ei, £2, £3, 
admits exactly four solutions. Indeed this equation expresses the condition of existence of a map of 
the 4-simplex in M'^ giving orientations e = ei, £2, £3 to the corresponding tetrahedra. Hence, first, 
it imposes G = 0, and therefore, according to the previous result /o4 = ^04- Now the data {/jj, /q^} 
determines the geometry of the complex in M^ modulo orientation reversal in this space (in other 
words simultaneous change of all orientations). Therefore for each ^q^, there exist exactly two sets 
of orientations r/e^ , r/ G {±1} for which the equation is satisfied, which proves the statement. 

Let (/o4!£i) denote one of the four solutions of G = u;=0. Then a variation (5/o4 of the length 
induces variations of 5G of the volume and 5u) of the deficit angle. Our third lemma, established 



in appendix 13 relates these variations to each other: 



dG 



doj\ 



04 

dG 



-2eie2e3 



V1V2V3 



04 



dl 



04 



T2V0V4/04 



(17) 
(18) 



The identity (|14|) is obtained by evaluating in two different ways the functional 5{G). We first 
consider G to be a function of /o4- We have: 



5{G)= J2 ^(l' 

eG{±l} 



04 



^04) 



dG 



dl 



04 



-1 



Eiirirvr'^(^o4 



eo.e4 



4/o4 V0V4 



'04 )• 



(19) 



In order to get the second equality we have used (jlSp . replaced the sum over the solutions by a 
sum over the orientations of the tetrahedra 0, 4, taking into account the overcounting by a factor 
2. 

Instead of considering G as a function of /o4; one can consider it to be a function which depends 
on the orientations ei, £2, £3 and the corresponding deficit angle uJq^. The following property 



G = 



3ei,e2,e3 such that Wqa = mod27r. 



yields: 



^iG) = \Y. '^('^04) 



(;i,e2,e3 



dG 



duj, 



04 



E 



ei,<:2,<:3 



l04 SJUJQ^) 
4 V1V2V3 



(20) 



The factor 1/2 in the first equality comes from the fact that, if uJq^ = 0, then co^^ = 0, so when we 
sum over all e, we are overcounting with a factor of 2. The equality between (J19p and (j20p leads 
to the identity p^ . 

It is worth noting that one can actually derive an identity slightly stronger that p4p. In fact, 
(jlSp divided by (|17p gives the derivative of the deficit angle with respect to the length: 



duj: 



04 



dl 



04 



, ,2 V0V4 



(21) 



It is then straightforward to check that 

/(^, 



Vr/, 



^04) 



V0V4 



1, ,2 '^(^04 - ^04) 



^^ol 



dl 



04 



/(/04) 



/± 



dZn4^i 



04*04 



6{uj\ 



rye- 
04 



V1V2V3 



fil 



04 J 



(22) 
(23) 



where the label appearing as an index of integration means that on integrates in the neighborhood 
of the value Iq^. 

In the next two parts of this section, we will take advantage of this identity to express the 
Feynman graph amplitude as a spin foam amplitude. 



10 



2. Spin foam in Feynman graph: inducing the geometry dynamically 

In order to understand the mechanism, let us begin by studying a simple example. Suppose 
that our Feynman graph F has the structure depicted in Fig[2| it admits five vertices and ten links. 
The two tetrahedra tq, t/^ form a triangulation of a 3-ball. According to the construction described 
in section |TI1 the amplitude of this graph is: 



foV(04) - - ° ^ 

where Or is the product of propagators 



feV(04) eo,e4 



fe)7^(04) 

The length of (04) depends explicitly on the nine other lengths. The previous identity (I14p allows 
us to remove this dependence by introducing a new label Zo4 and summing over it. Plugging it into 
([21]), and expanding the periodic delta function, as a sum 2t:5{uj) = Ylsez ^"^^^ ■, we get: 

^r = ^/ n h.dkj dlo4l,J2 E TrTn7^r(/.„/o4) (25) 

Notice that the tetrahedra Ti,r2,r3 form a new triangulation of the ball. In this expression, one 
begins to recognize the structure of the model (J14p . since the total action can be expressed in terms 
of the deficit angles associated with the edges of A 

S'a = E ^^"^^ = E '"^^^e- (26) 

T e 

Let us summarize where we stand so far in addressing the first step, which was to transfer the 
information about flat geometry from the measure into relations ll,{le)', it has led us from ([T]) to 
(jlOp . Now considering two tetrahedra sharing a face in the triangulation Aj., we have learned from 
the example above that the 'missing' distance variable can be replaced by a free label that we 
integrate over. The price to pay is the apparition of a constraint term - the delta function of the 
deficit angle - which, in the language of Regge calculus, imposes locally, the geometry to be flat. 
Ultimately, by expanding the delta function, we begin to see a model emerging from the Feynman 
graph, where the flat geometry is induced dynamically. 

The result (llip will be proved in two steps. First, we will to interpret our key identity as 
expressing the invariance under Pachner moves, and thus establish the topological invariance of 
the model ()12p . Secondly we will write the measure of the Feynman graph as the partition function 
of this model, computed on a specific triangulation. 

C. Topological invariance 

In this section we show that the spin foam model (I12p . 



11 



is independent of the choice of the triangulation and depends only on the topology of the manifold 
that A triangulates. The label (GF) appearing as an index of integration means that the integrals 
have to be gauge fixed, in order to properly define the partition function. 

The necessity of gauge fixing the model is manifest if one carefully looks at the integrand of (J12p . 
In fact, the sum over integers s of the exponentials leads to a product over edges of delta functions 
S{ijl). The constraints Wg = mod 27r are not all independent however, and the redundancy of 
delta functions induces divergences. These naive divergences are related to the existence of gauge 
symmetries in the model, which have to be gauge fixed. This is analogous to what happens in 
the study of the 3d quantum gravity amplitudes, as explained in [13). The symmetries and the 
corresponding gauge fixing of the model will be studied in great detail in the next section . 

In order to specify this gauge fixing, we first choose 4 vertices which belong to a tetrahedron 
To of A, and then assign to any other internal vertex u of A a tetrahedron r^, to which this vertex 
belongs. This assignment is said to be admissible if the following condition is met: given two 
vertices vi and V2 of A that are not in tq, then (^1^2) '^ t^i n t^^, that is, the edge (^1^2) (when 
it exists) does not belong to the intersection of the corresponding tetrahedra. We denote such an 
admissible gauge fixing assignment (tq, {tv}v^to) t>y T. We say that a tetrahedron r belongs to T 
if either r = tq or r = r^ for some vertex not in tq; we say that an edge e belongs to T if either 
e G To or e E T„, and if e admits v as one of its vertex - for each t„ there are three such edges, 
carrying the labels Ivi,lv2,lv3- The admissibility condition insures that these triplets of edges are 
independent, that is, we cannot have Im = ly/j. We will see later that the admissibility requirement 
can always be met. 

We now define a gauge fixing term 6qp and a Faddeev-Popov determinant Dpp as 

simile, Se) = ll6{le-l°){27r)6s^,so (28) 

±, ^ \l'vll'v2l'v3 / 

TfcJ 

The gauge fixing term imposes the values of le and Se to be l", s°, which are arbitrary fixed. Fixing 
the labels s's living on the edges of the assignment removes the redundant delta functions of the 
deficit angle at these edges. Now note that the constraints loI{1) = [2tt] act both on the labels 
Vs and the orientations. Therefore by removing them, we not only remove redundant constraints 
on the Vs but also necessary restrictions acting on orientations. As shown in appendix, one can 
reintroduce these restrictions by inserting an additional gauge fixing factor^ in Sqp'. 

The function x(^) is a characteristic function defined to be constant, with value 1, on 47rZ, and 
elsewhere; and il^ = '^^-^y ^t-^t is the algebraic sum, over all tetrahedra sharing the vertex u, of 
the solid angle A^ at v within r. 

The gauge fixed partition function is obtained by inserting in the integral (|27p the product 

"GF-^FP- 

Showing the topological invariance of the model amounts to showing the invariance of Z/\ under 
the Pachner moves (2, 3) and (1, 4), and its independence of the choice of T. The invariance under 



^ This factor acts 'on-shell' as Kronecker symbols for the orientations and is understood here as a gauge fixing. An 
interesting open question is whether this constraint can be expressed off-shell and implemented dynamically. 
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the move (2, 3) follows from the pentagonal identity, 

. , pJeoSo pie4,S4 1 ^ /" ^ p«<:i5i „ie2S2 aif^iSs 



<:0j«4 ei,<:2,<:3 



where tq, • • • r4 labels the five tetrahedra which triangulate the boundary of a 4-simplex, tj being 
the tetrahedron where the vertex i is omitted; the invariance under the move (1,4) is due to the 
(1,4) identity 

E V = (^ ^ / n d^o.^^, E n -^^o.D,, (31) 

(5gf denotes the gauge fixing term and Dfp the Faddeev-Popov determinant: 

<5gf = (2^)=^ n -^(^0* - ^o°i)'5so„sg,x(^^), D^y = 3' (32) 

where /g^, Sqj are any fixed values. At this stage the insertion of these terms are needed for topo- 
logical invariance; they will be fully and independently justified in the next section (together with 
the gauge fixing denomination). 



The proof of (j30p is a simple adaptation of the argument leading to the key identity (J14p . As 
before, we denote by r/e^ ,r/ = ±1, the two sets of orientations for which the deficit angle a;o4(/o4) 
vanishes. If this dihedral angle is zero modulo 27r, the 4-simplex [01234] can be mapped in M^. 
Each edge {ab) belongs then to a complex of three tetrahedra embedded in M'^. One can therefore 
find orientations Eq ,e^ of tq,T4, for which all deficit angles co^ are also zero modulo 27r. The 
actions of the (2, 3) move evaluated for Iq^ = Iq^ are then related by 

e^So + efS^ + efSt + efS"2+efSi^ ^ Sabu^f, = [27r] (33) 

{ab)^{OA) 

The superscript o means that so4 = in the actions. This equality, together with (p3|) . show that 

Vry, \ \ =/,dWoV(c^of)nh^ ^^^ 

Vo V4 Ji± fj^ H 

Now we want to sum over r/ = ±1. Since t/Ej are the only possible choice of orientations for which 
a;Q4 = admits a solution around Iq^, one can replace, in the RHS, the summation over 7/ by a sum 
over all orientations. Then summing the contributions of the two values l^^ leads to 






Now it is easy to check^ that the summation in the LHS can be replaced by a sum over the values 
of orientations £0,64. If one now expands the delta function in the RHS as a sum over so4 of 
l/27rexp (^8041^04)) one obtains (j30p . It's worth noting that we could have inserted any function 
/(/04) of the length /o4 in both sides of all our equalities, without modification, and therefore (f30ll 



In fact, our definitions of orientations (see appendix) imply that Eq ^4 = ±1- It shows that the application 
(±,7j) 1— > (t^Eq ,7764 ) is surjective onto { — 1, +!}■ 
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should be understood as an equality of measures on the space of such functions. Since ^4 and V4 
do not depend on /o4i we multiply both sides of ([M|) by V4e~^^'^4, ^i^ to get 

Vr?, '—^ = I d/o4 llA^T, ) n ' . % V|, (36) 

where in the RHS the values of l/a = l^^ are fixed to be any value - since the reasoning leading to 
(j34p does not depend on it. We now insert the trivial 'gauge fixing' identities 

1 = I \{dkdli{'-^^^^^ (37) 

1 = E7iEn'^^-.o°.(2-)'x(f^^), (38) 

where e = {rje^ , €4). As can be easily seen geometrically, and as more precisely shown in appendix, 
the algebraic solid angle fig, evaluated at Iq^, vanishes (modulo 47r) only if €4 = r]e^ , and therefore 
its image by the function x £^cts as a Kronecker symbol for £4. Hence, we obtain 



2 Vo (27r)4 



nd/oi/g,y^^d/o4/o'4E'^K'I )n^:^ ^.JgfI^fp. (39) 



for all r] = ±1. Eventually, summing as before both sides of this equality over 77, adding the 
contributions of l^^ and Iq^, and expanding the delta function in the RHS, show (|3ip . 

D. Symmetries and gauge fixing 

We have seen in the previous part that, in order to have a model which is topologically invariant, 
we need to insert gauge fixing and Faddeev-Popov terms in the definition of the spin foam model 
(J12p . We now want to show that these factors are not only needed for topological invariance, but 
arise naturally from the gauge fixing of the symmetries of the action. 

The action that arises from our construction reads: 

SA[SeJe] ='^erSr{Se,le) =^SeWe(^e) (40) 

T e 

The naive divergences of the model (J12p can be seen as arising from symmetries of the model 
mapping solutions to solutions. Fixing these symmetries, treated as gauge symmetries, will lead to 
the Faddeev-Popov determinants described in the previous section. In order to characterize these 
symmetries infinitesimally, we will study the zero modes of the Hessian of Sa- Here Sg will be 
treated as a continuous variable, in order to simplify the analysis^. 

Let us first describe the classical solutions. The equations of motion (eom) are 

= ^ = u;lil) Ve (41) 



Ssp, 




Ve (42) 



^ We can treat s^ as a continuous variable if one authorizes oje to be shifted by 2mr for n an arbitrary integer. This 
is clear from the expression of the periodic delta function S{uj) = ^n J dse^"^"~ "'^'. This affects the classical 
solutions but not the equations that involve derivatives of cu, which are our main interest here. 
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The first equation expresses the flatness condition. A set of lengths {l^} is solution if the complex 
can be locally* mapped in M^ in a way that respects orientations e. The second equation is 
reminiscent of the Schldfli identity for flat tetrahedra 

e'Gr e' 

This shows that a solution of (j42p is given by the values Sg = ale, with a an arbitrary constant. 
Notice that if one inserts this solution into the action (j4Up , one recovers the Regge action of discrete 
3-D gravity 

5^ = a5^/ew^(/e) (44) 

e 

In the following a solution {1^,3^ = l^} of the com will be called a Regge solution (we put the 
overall scale factor a = 1 for simplicity). 

Zeros modes are infinitesimal transformations of the fields (here the labels) that belong to the 
kernel of the Hessian 5'^S, computed on shell. This kernel is characterized by deformations 6se, 5le 
that satisfy the equations 

E|;«e. =OVe (45) 

e' 

E-"£e«e'+Et^*V^O Ve (46) 

e',e" e' 

where the label e has been dropped for clarity. This system of equations can be rearranged if one 
uses the derivative of the Schlafli identity 



^ 
5L 






SuJe . sr^ . b'^LOpii 



+ E'-e^ = ». OT 



(5L' ^—^ blebl 



and the fact that the matrix ( ■^\ is symmetric, as the Hessian of the Regge fonction (j44p . We 
can write (1460 as 

E (^^" - ^^") irS^ '^^^' + E ^ ^(^^' - ^<=') = 0- (48) 

e ,e" e 

If one focuses on fluctuations around Regge solutions, the system becomes equivalent to 

e' e' 

One can easily find infinitesimal variations ble-, 6se solutions of (j49p . Indeed, an embedding in M^ 
of the complex A, with the lengths {/g}, gives positions to the vertices. Infinitesimal moves of the 
vertices in M^ around these positions induces deformations of the edge lengths. By construction, 
the new deficit angles still vanish. Thus, these deformations correspond to one of the symmetries 



A global map exists, and is unique modulo translations and rotations, if the manifold that A triangulates is simply 
connected [l7| ]. 
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we are considering, acting on vertices of the triangulation. The positions of the vertices in an 
embedding are therefore pure gauge in our model. The generators are 3-vectors a^ attached to 
each vertex v which induce displacements of v in M^. They determine variations 5se = 6le{av, l^), 
51 e = 5le{l3v,le), which are solutions to the system of equations ([19|) . 

These transformations map solutions to solutions^; the corresponding symmetry group is 
(M^ X IR^)I^I, where |t;| is the number of vertices of the triangulation. The gauge fixing is per- 
formed by fixing a subset of edge lengths while taking into account the Jacobian or Faddeev- Popov 
determinant. These Jacobians can easily be read out from the results of section |II1 We start by 
choosing four vertices that form a tetrahedron tq in the triangulation. For these vertices we have 

d'a^d^a2d^asd'a, = 2d'adA^^^j^^. (50) 

The factor 2 is due to the sum over values of the orientation for tq. For each additional vertex v, 
we choose a tetrahedron r„, and denote l^^ , lv2 , U^. the three edges of r^ touching v. In this case we 
have 

d3a, = 2niil!^. (51) 

Each time we add a new vertex u„, one should insure that r^,^ is picked in such a way that the 
three corresponding edges do not overlap with any of the previously chosen edges l^^^^ ■ ■ ■ Iv^-u-: ^^ 
order for an edge to be selected at most once in this process. We are going to see in the following 
how this can always be concretely implemented for any closed triangulation. 

This procedure specifies an admissible gauge fixing assignment T^, as defined in the previous 
section, which allows us to gauge fix the symmetry acting on Vs. The gauge fixing term is given by 

(52) 
and the determinant by 

^^^^^ -d3«dAA.,^,d/e-' n.er,V/ ^''^ 

where \v\ is the number of vertices of A. The same procedure is applied to the gauge fixing of the 
symmetry acting on s's, with a gauge fixing assignment T^, and the corresponding gauge fixing 
term and determinant 

^GF= J{{27^)6s.,si\{x{^l{le)) (54) 

e&Ts veTs 

^"^^P^ -d3adAAeeT=d.e"n.eT^V/ ^^^^ 

There is no factor 2 in this Faddeev- Popov determinant, since here orientations of the s's tetrahedra 
are fixed thanks to the factors x(^) iii the gauge fixing term. The insertion of the functions x 



^GF = n ^(^- - 


-O 




Av d^«« 


= 2^"^- 


-3 HeeTi ^e 


d3adAAegT;d/e 


UreT, V. 



^ These symmetries, acting on the space of solutions, do not however leave the full action invariant away from 
classical configurations. Whether or not these symmetries admit an off'-shell extension is a interesting open question 
not addressed here. 
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restricts the values of orientations to those for which solutions of the eoni (j4ip , around which the 
zeros modes which have to be fixed, exist. One thus recovers the gauge fixing described in the 
previous section when one chooses the same gauge fixing assignment Ti = Tg for / and s. 

Finally, we need to show that an admissible gauge fixing assignment can always be reached. 
One way to construct such an assignment is to choose a rooted tree (also denoted T for simplicity) 
in the one skeleton of the dual triangulation. The one skeleton of the dual triangulation is a 4- 
valent graph; each vertex of this graph corresponds to a tetrahedron of A. A tree is a connected 
subgraph of this skeleton which contains no loop. A root is a choice of one particular end point of 
this tree diagram. We use the notation v* for the vertices of the dual triangulation, and r^* for the 
tetrahedron of A dual to i;*. A tree T is said maximal if each vertex of A belongs to a tetrahedron 
dual to a vertex of T: Vw G A, 3t>* £ T,v C t^* 

Given a maximal rooted tree T one can construct an admissible assignment as follows. The 
root of T corresponds to a choice of an initial tetrahedron tq. Then, lets pick an edge e* that goes 
from the root to a nearby tetrahedron Tg*: this edge is dual to a face /g* of the triangulation. If 
we consider Tg* — fe* , this consists of three edges meeting at a vertex Ve* G Te* ; in this way we 
define a correspondence between e* and a pair (re*,Ve*). This correspondence can be extended to 
the all tree if one uses the natural orientation of the rooted tree. Indeed, since the tree is rooted, 
every edge of the tree can be oriented with the orientation always going farther away from the 
root. And to any edge e* of T we assign a tetrahedron Tg* corresponding to the end point of this 
edge and a vertex Vg* being the unique vertex in Tg* — fe*- This assignment is clearly admissible. 
In order to see why, lets suppose that there exists an edge of A that belongs to both t^* — fe* 
and Te'* — fe'*- This edge would belong to the intersection; but this intersection is empty since 
tetrahedra always intersect on faces dual to the rooted tree. Eventually, the maximality of the tree 
implies the maximality of the gauge fixing assignment. 

E. Observables and partial gauge fixing 

In this section we will show our main result, namely that the original Feynman integral (jlOp 
can be expressed as an expectation value of an observable in our topological spin foam model. 

We have seen that the pure model ([T2|) has to be gauge fixed; the complete gauge fixing is 
performed by choosing a admissible assignment T, and by inserting gauge fixing and Faddeev- 
Popov terms K^^''^'{le,Se) = 5Qp5QpDppDpp^ with Tg = Ti = T. We will call observable of the 
model a function f{le, Se) of the labels, and define its evaluation as 

<«- = (^L n "''- ^ "-'""'< n — ij— pe) 

/ is not in general gauge invariant and therefore its insertion modifies the discussion of the gauge 
fixing. We will be interested in particular in the class of observables that are functions of the 
lengths only, of a subgraph F drawn in A. The typical example is the so-called matter observable 
Orile) = HeGr G^i^e), where F is a subgraph of the 1-skeleton of A. 

The introduction of such observables breaks a part of the gauge symmetries - namely the 
symmetry of the field le which acts on the vertices of F, turning the former gauge degrees of 
freedom into dynamical degrees of freedom. In the presence of matter observables the gauge fixing 
can no longer be symmetrical, and one needs to choose different assignments for the Vs and s's 
gauge fixing. Let Ts be a maximal admissible assignment {r^}, and Ti C Tg a subset that is a 
maximal assignment in A\F. In other words, Ti = {r^}t,^r- We then take care of the remaining 
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gauge symmetries by inserting the gauge fixing and Faddeev-Popov terms corresponding to these 
assignments. The evaluation of the observable Or in the (partially) gauge fixed model reads then: 



(Or>A = ^^^y n d^^^' E ^r(/e) f n —^^ ) K^"^'ile,s,) (57) 



We now suppose that F is a Feynman graph and A a triangulation of the 3-sphere S^ such that 
F C A. We want to show that the expectation value (j57|) is equal to the Feynman amplitude pO|) 
of the graph F. 

Recall that the triangulation A^ involved in (|lUp is a triangulation of a 3-ball, without any 
internal vertex nor internal edge; it is not a triangulation of S^. It is easy however to construct a 
triangulation of S^ from A^, simply by taking its double denoted -DA^ = Afc(J52Afc, built with a 
copy of Afc in the interior of a 3-ball B3 C S^ glued to a copy of A^ in the exterior 5^\i?3. 

The proof proceeds in three steps. Primarily, since both DA^ and A are triangulations of S^ 
which contain vertices of the graph F, the two triangulations can be constructed out of each other 
by a sequence of Pachner moves that do not remove the vertices of F. Secondly, (Or) A is invariant 
under these moves, provided that the variable le of each edge e G F that is removed by a (2, 3) move 
is replaced by its value 'on shell', namely the one fixed by the constraint ujI = 0. The justification 
comes directly from the analysis of section IIIICI One can therefore write 

(Or>A = (Or)DA,, with dr = Or(le,ll'{le)). (58) 

The quantities lli{le) are by construction the Euclidean distances, in any embedding^*^ of -DA^ in 
M^, between the vertices of F that are not connected by the edges of the triangulation. 

Thirdly, /r = {Or)DAk turns out to coincide with the Feynman amplitude Ir given by ([TO]) . 
In order to see this, let us express the gauge fixing terms for the triangulation DAk- We need to 
choose a maximal assignment Tg of DAk and a assignment Ti C T^ maximal in DAk \ F. Notice 
first that, since all the vertices of DAk belong to the graph F, Ti is empty. Next, recall that, as 
emphasized in part [Til the triangulation Ak of the ball is dual to a tree T; we can take this tree 
to be a gauge fixing tree Tg, according to the procedure described in IIIIDI Vertices, edges and 
simplices of Tg are the vertices, edges and simplices of A^ C DAk- Therefore the gauge fixing is 
performed by inserting 

€'f = n (2vr)5s.,sg n ^(^^) ; ^FP = TT^ ^^^) 

eGAfc i>gAfe -l-leeAfe e 

Since Ak has no internal vertex or edge, A^ and DAk possess the same number of vertices and 
edges, all located on the boundary of the ball B, whereas the number of tetrahedra in DAk is 2k. 
For each tetrahedron r in the interior of B, there exists a copy of r in the exterior of B, sharing 
its edges with r, and consequently having the same volume. Moreover, the tetrahedron r and its 
copy have independent orientations e,- and e'^. Therefore, using (f59l) . /r can written as a quantity 
computed on the triangulation A^ 




Ir = y n ^^^^^ E ^r(^e,o' (le)) I n ^ I e^^^^^^-^ n xi^r) m 

eSAfe e,e'e{±l}*^ 



^ Note that DAk can always be embedded in R . 
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where deficit angles are given by 

^r'ile)=Y.('r + e'M (61) 

Finally it's easy to convince oneself that the total algebraic solid angles seen from vertices v vanish 
modulo 47r only if e,- + e!^ = for every tetrahedron of A^. Therefore Jp = -^r- 

IV. ALGEBRAIC STRUCTURE 

In order to complete our investigation, we now want to give an algebraic interpretation of the 
topological state sum model introduced before. This will show that this model is really a spin foam 
model, namely constructed purely in terms of algebraic data. 

A. A Poincare model 

1. The result 

We first introduce some notations. We define 

p={le,Se), J dfi{pe) = ^^ J dljl (62) 

and introduce a symbol associated to the tetrahedron r with vertices 0, 1, 2, 3: 

P23 P13 P12 1 ^ V- e"-^-(^^.^''^-^) ^ COS Sr {S I J, llj) 

Poi P02 Pm\~ ^ V{lij) V{lij) ' ^ ^ 

with 

Sr{le,Se) = Y,^e9l{le)- (64) 

eSr 

The previous model can be written as a Ponzano-Regge like model |l9l | 

Za= [ll^piPe)Ul''^' ''^' ''^^l' (65) 

7 V r I ^^I P'l ^4 J 

where the integral is over all edge labels and where ej label the 6 edges of r. 

The topological invariance of this model is essentially due to the identity (|30p . which can 
rewritten in terms of this symbol as a pentagonal or Biedenharn-Elliot identity 

PI P2 P3l fpi P2 P,\ U„^p)[P^ P^ P^\[P^ P^ P^\[P^ P^ P^\. (66) 

P4 P5 /56 J I P4 Ps Pe J J I Pe P P5 J I P4 P Pe J I P4 P P5 J 

All these considerations clearly suggest an interpretation of this symbol as a 6j-symbol for a 
classical group. This is indeed the case, according to the following 

Theorem: p = {l,s) labels the unitary representation of the 3- dimensional (Euclidean) 
Poincare group IS0(3) with mass I, and spin s; dp{p) is the Plancherel measure which arises in 
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the decomposition o/L^(ISO(3)), and the symbol {pi} defined in I163\) is the normalized Qj -symbol 
o/ISO(3). 

As a first corollary of this theorem, we get an other independent and non geometric proof of the 
pentagon identity (j30p . since this is a direct consequence of recoupling theory. An other corollary 
of this theorem is the fact that this 6j-symbol satisfies the orthogonality relation 

mpu) I ^- ':' ':' ] I ^- ':' 7 ] = 2.6'-^^. m 



POl PQ2 P03 J 1 POl P02 Po3 J ''"^''*"^ ^03 

It is interesting to give a direct and independent geometric proof of this identity. We consider two 
tetrahedra r, r' with vertices 0, 1, 2, 3, with edge lengths and spin labels {hj, sjj), {I'jj, s'jj), these 
labels being such that Ijj = I'jj, sjj = s'jj except for (IJ) = (03). We want to compute 



(68) 



^ = ^2^7 ^^12^12— ^r ^j;— 

= fdlulhT^^-^^^fv^^i n e^(""^^--'<^OJ (69) 

•^ e,e' ^ ^' \(a6)^(12) / 

where Oij denotes the dihedral angle at the edge (ij). One can now use (|A14p 

to express the change of integration variables 

dlulu = d9o3^ (71) 

'03 



and 

9^12 



^^(£^12 - e'e[,) = <5,,,,5(/o3 - ^03) 



dl 



03 



Se,e'S{l03-l'03)j^ (72) 

'12'03 



Since Zo3 = ^03' ^^^ ^^^ volume factors in the integral cancel out; and since Sab = s^j, for (ab) ^ (03), 
the integral reduces to 

/ = ^i^^ r d^03 E e-^-(--^03) = ^(^03_i^2vr5,,3,,3 (^3) 

'03 ^0 ^ '03 



which is the desired result. 



2. Proof of the Theorem 



We will need some definitions. We start by describing the representation spaces of the Poincare 
group. Let us fix a frame {ex,ey,ez) of R'^. The carrier space corresponding to a representation^^ 
p = (/, s) of 150(3) = SO{3) X M3 is given by: 

£p = We C\SU{2))/y^ G [-7r,7r], cl){xh^) = e^^«</.(<7)}, (74) 



The spin s € |Z label representations of the little group (7(1). In the following we restrict to the case where all 
spins are integers. 
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where h^ = e*^'^^ " are elements of the [/(l)-subgroup of rotations with axis z and angle .^. The 
action of a rotation g G SU{2) and a translation a G M'^ on this space is then 

{g, 0) ct){x) = (P{g-'x), (1, a) 0(x) = e^^-'^ (/.(x) (75) 

where P^. = a; • / is the image, by the rotation associated to x, of a fixed vector I = le^ of length I. 
We now define the Clebsh-Gordan coefficients. The tensor product (/i, si)® (I2, S2) decomposes 
on representations p^ = {l^-,s^) such that li,l2,h satisfy triangular inequalities and such that the 
sum of Sj's is an integer 

\h - hi < h < h + h, si + S2 + S3 eZ 

Given three lengths satisfying these inequalities, one can construct an oriented triangle in the (x, z) 
plane such that 

^3 = aei3 • h + 0023 ■ ^2 (76) 

where one introduces the 2x2 matrix representing a rotation of angle and axis y 

_ f cos 61/2 sin6'/2 
"^~ \^-sin6'/2 cos 0/2 

The angles ^13,^23 G [0)'^] are the angles between the edges (1) and (3) (resp. (2), (3)) within the 
oriented triangle [123] whose edges lengths are li,l2,h- They are characterized by 

/2=/2_^/2_ 2^2/3 cos 6*23 (77) 

and similarly for ^13. The Clebsh-Gordan coefficient: 







also called 3j-symbol, is defined (up to a normalization factor) to be the invariant map Si^^si ® 
£i^s -^ ^ik,Sk- "^^^ structure of this object is similar to the Clebsh-Gordan coefficient of the 
quantum double VSU{2)^ whose construction is described in [20[]. When properly normalized, it 
can be written as 

— 1 3 

CS^^ = (^) ' /^^^^ n(dC.e-'«0 '^.,.-(^3/^C3«^»)'^.,.-(^3/^C3%;) (78) 

where dQ is the normalized Haar measure on [/(I). It is clear that this coefficient is such that 
<^rf ^t h . h = C7^ig2S|e-<"i^i+"2'^2-s3C3). Moreover, it is straightforward to check (thanks to 

(|76l) ) that it defines an invariant tensor, that is, Cgx^ = C§1, and that X^j Xj • li=0. The condition 
that imposes the sum of the spins to be an integer comes from the identities Cx] = C^^. = C^'^ = 
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The normalization factor is chosen in order to have the resolution of the identity [3] : 

I d^(p3) (/ dx^CPl^if^Cil^^ = 611 i^iKl (^2) (79) 

where the 'bar' denotes the complex conjugate, dfi^p^) denotes the Plancherel measure ([62]) . and 

Ju{i) 

It is worth mentioning two further properties that the Clebsh-Gordan ()78p satisfies. The first 
one is a 'braiding' property which describe the change of the Clebsh-Gordan under the interchange 
of the labels i and j: 

c£r.\z = i--^r^''-''cpi^r.i^ (so) 

derived thanks to the identity /iTrOe^-Tr = O-e- The second property expresses the dual intertwiner 
CP1P2P3 : s^^ ^ s^^ s^^ in terms of the 3j-symbol ([78]) : 



Finally, it will be important to note that, if one restricts, as we do, the spin to be integers, one 
can restrict the integral in ([78|) to be over ^ G [0,27r], provided that one takes the SO{3) delta 
function in (ffSf) . namely Sg{x) = ^{6g{x) + 6-g{x)). 

One can now construct the 6j-symbol, defined as a contraction of Clebsh-Gordan coefficients 

/ rlxmrl'rnnrl'rno nP^^PO^POl (^P03pi3p0l rip02p23p03 = J ^23 Pl3 Pl2 I ^P23P13P12 ,'oo\ 

/ aXoiaXo2aXo3(-^a;i2a:o22;oi^xo3a;i3a;oi'-^a;o23::23a;o3 — I p p p f '-"X23a;i3xi2 • V°^J 

The way the contractions occur is dictated by the geometry of a tetrahedron (0,1,2,3) with ori- 
ented edges (/J). There is one Clebsh-Gordan attached to each face of the tetrahedron, and one 
representation for each edge. 

Let us now prove the theorem, by first noting that each Clebsh-Gordan of the LHS of ()82p 
supplies two (5-functions on SU{2), so the integral contains 6 delta functions. The integration 
is over three SU{2) elements xoi, and 6 U{1) elements, denoted by ^oij^oi- -^y performing a 
simple change of variables XQi — > xoi/i^g. , one sees that the integrand depends only on the xqj 
and the differences ^oj ~ S,'oi- Then the integration over the SU{2) variables allows us to solve 
straightforwardly three of the delta functions. In order to write the result of integration in a 
convenient form, we assign an angle ^/j to each edge (U), < /, J < 3 of the tetrahedra [0123], 
and denote the corresponding group element by hjj = hjj = h^jj. We also denote by a/^ = aj^'^ 
the rotation clq associated to the geometric angle between the edges (/J) and (IK) in the 

triangle [UK]. The later angles are fully determined by the lengths Ijj. After integration, the 
LHS of (1H2D takes the form 



n [d^ije-''"^") S{hoia^^ho2a^^ho3a^')5^^^^-i{xuhi2gi)S^^.^^,-i{xuhi2g2) (83) 



KJ 



times the prefactor (vr ^^01^02^03) ^(^r ^^12^23^31) ^ , where we have defined 

g, = {a?^)-'h,oa^'hos{a!^Y\ 92 = (afO)-i/iio(af )-i. (84) 



22 



One can simplify the expression for gi if one remembers that the sum of the three interior angles 
of a flat triangle is equal to vr. Using this fact for the triangles [012] and [023], we get (a{'^)~"^ = 
a-T^a2^aQ^ and (03^)""'^ = a_^ao^a2^. Plugging this into the expression for gi yields 

gi = a-T^a2 {a^ /ioi«o ^os^o ^02} ^20 '^2 0-71 = 0^7^02 /i2o ^2 o_7r = — 02 /i20O2 , (85) 

where for the last two equalities, we first have used that the term inside the bracket is 1 - this is 
due to the first delta function in ([83]) -. and then the property a-j^ha-,^ = —h~^. 

In order to complete the proof of the theorem, we first have to find the solutions of the constraint 
expressed by the first 5- function in (jSSp . These solutions of the constraint arise from the flatness 
of the tetrahedron, which leads to a remarkable identity at each vertex (see (IBlh ): 

h\jai^'h\Lat''h\KafJ = -el V/ (86) 

{IJKL) is an even permutation of (0123), Kjj = h^e is the rotation associated to the exterior 
dihedral angle of the edge (U). Furthermore, if one considers (|86p as an equation in S,ij, then 
CiJ = ^fj) and this is the unique solution belonging to the interval [0, 27r]. Multiplying ([86l) by 
Ctt on the left, and a_7r on the right, has the effect to change all h into h^^, while leaving the a's 
unmodified; is one now multiplies the LHS by /i2^ = — 1, one sees that ^jj = 2tt — (p'fj is a solution 
of (|86p with +1 in the RHS. Theses identities are carefully studied in appendix [Bj The bottom 
line is that the equation /ioiod^^020o^^030o^ = if (recall that we work with integer spin) admits 
two solutions when ^jj S [0,27r], namely 

Cfj = <Ph, or Cjj = 2vr - 0f^, (87) 

We denote the respective h^± by hjj. 

For each of these solutions, we use ([86]) to express the Euler decomposition of gi and (72 

g^ = {a^^r'hUa?')-' = -e{h\2)-\!\h\.^-' (88) 

92 = -a^"hl,a^^ = eihl,)-\at)-\hl,)-' (89) 

Using this decomposition one can rescale ^ij — > S,ij — ^.L, I < i,j < S in ([83|) . and factorize from it 
the Clebsch-gordan coefficient Cxllxllx^i^2 ■ '^^^ proportionality coefficient for each set of solutions 
is 

n '^~"'"'^" ( W^) / ( n d^o*^"^"' ) '^(^oia(}'/io2ao'3;^o3ao'^) (90) 

where the subscript e in the integral reminds that we need to integrate around the solution ^q? = Coi ■ 
Now if we denote g{Coi) = /ioiad^/i02ao'^/io3) (|B14p gives the Jacobian of transformation from 
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the Haar measure to the S,oi measure around a solution 5'(Coi)'^o^ = =^1: 



dg = 7rp^lld^0^ (91) 

'01t02t03 

I 

where the volume of the tetrahedron appears. Therefore the expression in ()90p can be written as 

n -"'«' idk) I (n ^««.--) ^««. - ««.)^ - %^^ (-) 

The 6j-symbol is eventually obtained by summing the two contributions e = it. Therefore it reads 

/ P23 Pi3 P12 1 _ . -,^y,ST.,n ^^^^r{sij,lij) . . 

I Poi P02 P03 / - ^ '^ ■ ' V{lu) ^""'^ 

The exterior dihedral angles have been replaced by the interior ones 9^ = vr — 0*^, leading to a factor 
(— l)^'*^-^ in the formula. 
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B. Poincare versus Ponzano-Regge 

We have shown in the previous sections that a Feynman graph amphtude is the evaluation of 
the observable Or{le) for a Poincare spin- foam model. On the other hand, we know [8| that a 
usual Feynman graph is also the abelian limit of the evaluation of the same observable for the 
Ponzano-Regge model. This means that one should be able to relate in a direct way our Poincare 
model with a suitable semi-classical limit of the Ponzano-Regge model. 

In order to do so, let {ue, Ve) be an edge-labelling of a triangulation A, with pairs of half-integers. 
We consider the square of the Ponzano-Regge model: 

^r'= E Il^uAv.ll^u}r{v}r (94) 

where {u}r and {v}r denotes 6j-symbols of SU{2) and dj = 2j + 1. After a simple change of 
variables 

u = j + s, V = j — s 

the partition function reads 

^r' = E 11(4 - (2^^)') U^J + ^}r{j - S}r (95) 

\Se\<je 

Next, let us introduce a length variable I = djlp, where Ip is the Planck length. We would like to 
evaluate the limit /p ^ of the model, in a sector where the typical value of / is kept fixed, and 
the spins s are bounded by a cut-off n much smaller than the typical value of j: n <^ l/lp- 

Interestingly, using the techniques presented in [l5|], asymptotics of 6j-symbols in (f95]l can easily 
be evaluated in this limit. We get: 



{±±S-i}2~l?-l-COs2 

^2lp 2^^ TT V{lij) 






(96) 



for fixed values of the length I and the spin s. This suggests the following semi-classical limit for 
the model: 

1 /73\ (I'^l^l^l) ^ 1 f JSr(lp) 

J- / S \ „Pr2 i / -TT ,, -rT,9 V^ V^ TT- e ^'•^> 



where Vr is the volume of the tetrahedron r with edge lengths /g; and the action is given by (we 
write 6*6 for 0g and e, e' for e,-, e^): 

^r%) = E ^(^ + e')^e0e(/e) + (e " e')Se0e(/e) + (e + ^)\ (98) 



eer ^^P 



We can reorganize the sum over the orientations 



12. 



{e^,e'^} V^=±l <=±1 



using 2 cos ^^^ cos ^-j^ = cos a + cos b. 
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and identify both Regge action and Poincare action. Denoting the RHS of (I97p by Z/\{n,lp), we 
are interested in the quantity 

Za = lim hm Z/\{n,lp) (100) 

We will show that Z/\ is the partition function of our Poincare model (J12p . Let us focus on the 
first limit Ip ^ 0, n being fixed. The integrand of Z/\{n,lp) is a product of sums which one can 
develop. The integral splits then into 21"^! terms of the form 

^ ' •' eGA {se<2n} r ^ rdKx \r]'^ ) T&K2 \ Vt / 

Ki and K2 are disjoint sets covering the all set of tetrahedra; \K2\ denotes the cardinal of K2. Note 
that we have made the change of variables Se — > 2se, and then the spin are now integers. All terms 
such that K2 7^ are oscillatory integrals, and therefore under stationary phase evaluation are 
typically of order Ip . The term corresponding to K2 = is independent of Ip and is therefore 
the only term surviving when Ip goes to zero. 

Eventually, by taking the limit n — > 00, we get the desired result and recover the Poincare model 
in this semiclassical limit. 



V. FEYNMAN DIAGRAMS ON HOMOGENEOUS SPACES 

In this section we consider a scalar quantum field theory on spherical space-time. A Feynman 
amplitude of a graph T embedded in the unit 3-sphere S^ takes the form 

Ir= f dui---duN n '^™(^ii) (103) 

•^^' (ii)er 

dui is the normalized measure on the 3-sphere. The integrand is a product of propagators, which 
are functions of the dimensionless spherical distances lij G [0,7r] between the vertices, and invariant 
under the action of the group 50(4). 

The results obtained for the flat case can directly be extended to the spherical one. Again, 
the idea is to gauge out the integration measure. The invariant measure is expressed in terms of 
spherical lengths and the Haar measure dU of 50(4) in the same way as in section HIl 

dui • • • du3+k = dU Y, n d^e sin /e n ^ (10^) 

eG{±l}''' eGAfc rGAfe '^ 

Afc is a triangulation, with spherical tetrahedra, of a 3-ball in S"^, which possesses 3 + fc vertices 
on its boundary, and which does not contain any internal edge. Vt is the square root of the Gram 
determinant det [cos kj] (see appendix) associated to the tetrahedron r. 

A key identity relates the Gram determinants of the tetrahedra of a spherical 4-simplex, which 
has a form analogous to ()14p . We refer the reader to the appendix where the formula (J17p and (jlSp 
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are established for spherical simplices. This leads to the emergence of the topological spin-foam 
model 



The action for each tetrahedron reads 



Sr = Y.Seei{le) (106) 



eSr 

where O^, is the spherical interior dihedral angle of the edge e in r. The Feynman amplitude /r is 
then the expectation value of the observable Or(^e) = HeGr G'm(^e) for the model (jlOSp computed 
on any triangulation A of S^ which contains F as a subgraph. Analogous results for hyperbolical 
space are obtained by working on the hyperboloid and by replacing all the angles by hyperbolic 
angles. 

Again, the model (jlOSp takes the form of a Ponzano-Regge model. The weight associated to 
this model satisfies a pentagonal identity, and can be interpreted as a 6j-symbol. It is not difficult 
to convince oneself that the underlying group is the quantum double T)SU{2) [20[|; the proof of 
this statement is identical to the flat case, so we do not repeat it here. 

The classical study of the action (|106p differs from the flat case, since the eom non longer admit 
the so called Regge solutions Se = ale, unless a = 0. Therefore the connection with usual Regge 
gravity seems to be lost. However, it is worth noting that these solutions can be reintroduced by 
adding a volume term 2^^ VolT-(/e) to the action (jl06p . where VolT-(/e) is the spherical volume of 
the tetrahedra (which should not be confused with Vr). 

One can check that the model obtained is still topological, and that Feynman amplitudes can 
be written as evaluations of observables for this model as well. Let us introduce length variables 
le = vALg, where A is the cosmological constant. The new action can be written in terms of the 
geometrical data of a tetrahedron of curvature A 

Sr, = Y, SeO:^ {Le) + 2K^'^\o\r^ (Le) (107) 

eGTA 

where t\ is the image of r by the re-scaling l^ ^ L^ = le/V^., and 91^ is the dihedral angle of e in 
T\. Now using the spherical Schlafli identity for the tetrahedron ta 

Y^ LedOe + 2AdVol = 0, 

e 

we see that the variation of the action (jl07p gives Regge solution Sg = vALg. If we insert this 
solution into the action (jlOTh . we recover the Regge action describing discrete 3d gravity with 
cosmological constant A. Thus, the presence of a cosmological constant removes the degeneracy of 
the Regge solutions, and leads to the natural value a = VA for the scale parameter. Eventually 
let us remark that, since the topological model associated to the action (jl07p reduces to the 
Poincare model (J14p when A ^ 0, we expect the associated symbol to be, again, the 6j-symbol of 
a deformation of the Poincare group. 

VI. GRAVITY AND BF THEORY 

Before concluding, we would like to give some continuum explanation for the results we have 
obtained here. The fact that our state-sum model is built from 6j-symbols of the Poincare group 
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shows, from standard arguments [1]|, that the topological theory which we have described in terms 
of a spin foam model is the quantization of a Poincare BF theory. This theory can be written in 
terms of a one-form E and a connection A, both valued in the Poincare algebra, whose generators 
Ji, Pi, i = 0, 1, 2 satisfy 

[J',J^]=2e*^'Vfc, [J\P^] = 2e'^'Pk, [P\P']=0, (108) 

and 

E = e*P, + b'Ji, A = c*P, + u'Ji. (109) 



The action is given by 



S= [ tr{EAF[A]) (110) 



where F{A) = dA + A/\ A, and where the invariant trace is taken to be tr{JiPj) = 5ij, tr(PjPj) 
tr(Jj Jj) = 0. In components this reads 



Jm 



^R\uJ]+b^^d^d). (111) 



This model was introduced in 21[ and recently studied in [23]. e is interpreted as the three 
dimensional frame field, uj as the spin connection; 6, c are some topological matter fields. Its 
equations of motion when 6 = c = have, as the only solution, the flat space-time d^e = and 
R{ijj) = 0. Following [^l; one also expects the Feynman diagram observables to be related to Wilson 
lines observables for this theory 



Or = Pexp / tr(eiJo). (112) 

This result is reminiscent to the Polyakov reformulation of Feynman amplitudes as worldline 
integrals [23] • Such an interpretation provided a natural way, by replacing worldlines by worldsheets 
integrals, to deform field theory in terms of a dimensionfull parameter - leading to string theory. 
Analogously, our formulation of Feynman amplitudes in terms of spin foam models may provide a 
new way to think about consistent dimensionfull deformation of field theory structure, in order to 
take into account the quantum nature of background geometry. 

VII. CONCLUSION 

We have given in this paper a new perspective on closed Feynman diagrams in three dimensions. 
The main point was to show that the language of spin foam models - which was developed in order 
to address the problem of a background independent approach to quantum gravity - naturally 
arises in the context of quantum field theory, and that a careful study of Feynman amplitudes can 
even lead to a purely algebraic understanding of the quantum weights. It was also shown that 
usual field theory can be given a background independent perspective in the sense that the flat 
space geometry is purely encoded in terms of a choice of weights controlling the dynamics of the 
geometry . 

This result is by itself very encouraging since it clearly shows that the language of spin foam is 
not unrelated to the usual language of field theory and also suggests a new way to investigate the 
semi-classical limit of quantum gravity. Namely, the Feynman diagram spin foam model provides 
a 'zero order' approximation of the full quantum gravity model, giving strong constraints on the 
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possible quantum gravity extensions. Such consistent extensions are known in 3d; but it should 
be clear from the analysis presented here and as shown in [l] that the results of this paper can be 
extended to higher dimensions. 

We have addressed in this paper the case of scalar field closed Feynman diagrams. This clearly 
needs to be investigated further; first to the case of open diagrams but also to the case of higher 
spin fields. We expect in the later extension, the variable s which appears in our model, to acquire 
a more direct physical meaning. It would also be useful to have a deeper understanding of how 
usual Feynamn diagrams can be understood as the expectation value of topological observables 
from the point of view of the continuum, as sketched in the last section. 
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APPENDIX A: GEOMETRY OF SIMPLICIAL COMPLEXES 

In this appendix we first recall some definitions and notations, and give a proof of the results 

1. Basic definitions and discussion about orientation 

We start by recalling basic features about D-dimensional simplices. An abstract (or combinato- 
rial) D-simplex is a D + 1-elements set {0, • • • , D}. Its k + 1-elements subsets are called /c-faces of 
the simplex. An embedded spherical D-simplex is defined hy D + 1 normalized vectors (eo, • • • eo) 



in M^+^, e? = 1. The spherical edge lengths lij are given by 






vjCo fjri -1 Ol * 0-7 • 






We denote by (bo, • • • , bo) the dual vectors: bi.ej = 5ij, and we define : 






V(eo,--- ,eD) = |det(eo,--- ,eD)\ 




(Al) 


G{eo,---eD) = det(ei.ej) = det(cos/ij) = V^(eo,- • 


• ,eD) 


(A2) 



V^ = G is the determinant of the Gram matrix associated to this simplex. The interior dihedral 
angles are denoted by 6ij £ [0, vr] and defined by 

bi.bj = —\bi\\bj\ cos 9ij 

In this paper we work only with the interior angles and refer to them simply as dihedral angles. 
9ij is the angle between the {D — l)-faces i and j (i.e the faces obtained by dropping the points 
i or j) ^^. Hence, they are labelled by the {D — 2)-face Fij opposite to i and j. In particular for 
D = 3, dihedral angles are labelled by the edges of the tetrahedron, and for Z? = 4 by the faces of 
the 4-simplex. 

An orientation of an abstract D-simplex is a choice of an ordering (0, • • • D) of its points, up 
to even permutations ^]. A coloring {le,€) of an oriented simplex is a set of numbers le G [0, tt]. 



The exterior dihedral angle given by (pij = tt — 6ij is the angle between the normals to the (D — l)-simplices i,j. 
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labelling the edges of the simplex and satisfying triangular inequalities, together with a sign e = ±1. 
A given coloring promotes an abstract oriented simplex to an embedded spherical simplex; namely, 
it defines, up to rotations in SO{D + 1), normalized vectors (eo,---e£)) - ordered according to 
the orientation - in R-^^^, such that coslij = Cj • ej, and such that the sign of the determinant 
det(eo, • • • en) is e. The embedding is then called compatible with the coloring. In the flat case, a 
coloring (/e, e), with le £ M+ and e = ±1, specifies the geometry of a flat simplex (0, • • • D) in M , 
namely a set of D vectors li, i = 1- ■ ■ D represented by the oriented edges (Oi), and such that the 
sign of the determinant det(/i, • • • Id) is e. For instance, 4 ordered points (0, 1, 2, 3), six positive 
numbers lij, i,j = • • • 3 and an element e G {il}, determine the geometry of a flat tetrahedron 
inM^. 

An orientation of a Z)-simplex induces an orientation for each of its {D — l)-faces by dropping 
the missing point in an even ordering of the points of the D-simplex in which this missing point is 
first. Now lets consider two oriented Z?-simplices sharing a (D — l)-face. The orientations of these 
L'-simplices are said consistent with each other if the shared face inherits opposite orientations from 
them. A triangulation A is said orientable if there exists a choice of consistent orientations for all 



simplices 2j]; such a choice is then, by definition, an orientation of A. In this work we consider 
only orientable triangulations. For each of the triangulations we are working with, an orientation 
is implicitly chosen and fixed; and we don't refer to it anymore. A coloring of a triangulation is 
then a set {le,ea}, where le G [0, vr] label the edges and e^ £ {±1} label the simplices. Abusing 
terminology, the signs e^ will be called 'orientations' of the simplices. 

With these definitions it is now straightforward to check the property illustrated in Fig[Tl section 
im Given two D-simplices ai and (72 glued along a {D — l)-face, a coloring Co-iUo-2 = {^ei £0-1, £0-2} 
defines a map 4'aiUa2 ~^ ^^ such that (pa^ is compatible with the coloring Ca-^ = {le,ea-i} induced 
on ai by Co-iUo-2- Then, if e^i = €^2, the two points opposite to the common face in each embedded 
L'-simplex do not belong to the same half-space bounded by the hyperplane spanned by this face. 

Let us give some precisions about the notion of deficit angle: we consider a set T,f = {a D F} 
of D-simplices sharing a given (D — 2)-face F in a colored triangulation. 6p denotes the dihedral 
angle of the face F in the simplex a £ T,p. The deficit angle of the face F is a function ujp{le) of 
edge labels and orientations e = {eo-,cr G T,p} of the D-simplices surrounding F. It is defined by 

^HQ = Yl ^-^F (A3) 

Deficit angles are usually defined modulo 27r. They represent the curvature localized on (D — 2)- 
faces, in the sense of Regge calculus. While each Z?-simplex can be mapped in S (or M in the 
flat case), in general there is no global map, compatible with the coloring, of the complex Sj?. 
The existence of such a map requires the deflcit angle ()A3P to vanish modulo 27r. We represent 
below, for D = 2, and in the case of three triangles surrounding a vertex, different configurations 
for which the deficit angle at this vertex vanishes modulo 2tt. The complex formed by the triangles 
is embedded in the Euclidean plane (in the flat case). We denote by e the triplet (ei, 62, €3), where 
ej is the orientation of the triangle opposite to the point j. 

We eventually define and describe some properties of the solid angles seen at vertices of a 
triangulation. In a 3L'-triangulation A, we consider a vertex v and a tetrahedron r to which v 
belongs. A 2d sphere surrounding v intersects r along a spherical triangle: the angles of this 
triangle are the dihedral angles 9^ of the edges of r meeting at v, and its area, denoted by At, 
is the solid angle seen at v within r. The spherical triangles associated to all tetrahedra sharing 
V triangulate a surface called the link L^ of the vertex v. The total algebraic solid angle at v is 
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£ = ±(1,1,1) 



e = ±(1,-1,1) 



£ = ±(1,-1,-1) 



defined to be 



Of, 



E^'r 



tDv 



1. Note that by construction Ql 



61 if e-r = 1, and 27r 



0,% mod 47r, 
% if er 



where ^^^ = yl,- if e,- = 1 and Att — Aj- if e^ - 
where the RHS is defined as in section [III CI 

We also define the quantity Col — Z^rDe ^e r; where 9, 
Note that this quantity is related to the deficit angle as tDg = Wg mod 27r. 

Using the relation between area and angles for a spherical triangle, which reads here Aj- 
^eDv ^e ~ ^5 it is straightforward to check the following equality: 



-1. 



1 

2^ 



o: 






■'ie) 



X{Lv) 



(A4) 



where x(-^) = I''"! ~ l-^l + |c| is the Euler characteristic of the link, |r|, \F\ and |e| being the number 
of tetrahedra, faces and edges of A touching the vertex v, or equivalently the number of triangles, 
edges and vertices of the triangulation of the surface L^. If A triangulates a manifold, every link 
Ly is homeomorphic to a 2d-sphere, and therefore xi^v) = 2. The equality (|A4p relates the total 
solid angle at v with the deficit angles of the edges meeting at v. 

Lets study the case of a vertex surrounded by four tetrahedra 1,2,3,4. This configuration 
arises, for example, after a move (1,4). First, if the deficit angle cJqj of the four edges (Oz) vanish 
modulo 2tt, then the Footnote HI together with (|A4p . yield: 



Q.Q = 47rn(e), with 



if 



n e 



y{i,j) 



if not 



(A5) 



In particular, f^Q = mod 47r. Now lets only suppose that the deficit angle oJq^, with e = (ei, £2, £3), 
of the edge (04), vanishes modulo 2it. With the arguments exposed in section [III C\ this means 
that the complex of four tetrahedra 1,2,3,4 can be mapped in M^, and therefore, there exists a 
value of the orientation el such that the three other deficit angles vanishes modulo 27r as well, 
which implies Oq" * = mod 47r. Now as a consequence , one has JIq" * 7^ mod Air. This 
justifies the identity (j38p : the insertion of the function xl^o) ^^^^ ^^ ^ Kronecker symbol (5e4,e°, 



where £4 



v4- 



Except for pathological configurations whose contribution is negligible. 
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2. On the geometry of the simplex 

We work here in the context of spherical geometry. The resuhs for flat space are obtained by 
taking the limit where lij -^ while ratios hj/hi are kept fixed. In this limit, G reduces to the 
square of the usual Euclidean volume, and we have the following correspondence : 



V(eo, • • • , cd) ~ D]V{lij)), sinlij ~ kj, cosikj) - 1 ~ --ifj- (A6) 



1, 





(A8) 


V(efc,--- ,60) 
Vibo,--- ,bk-i) 


(A9) 


Vn) within span{vi, ■ ■ 


• ,Vn). 



Abusing terminology, we will call V the 'volume' of the simplex, keeping in mind that it reduces 
to the volume only in the flat space limit. 

Next we will need the following result. Let (eo, . . . , ejj) a D-simplex and L^ be the subspace 
span{6o, • • • , &A;-i} generated by the vectors dual to (cq, . . . , e,fc_i). By construction L^ is the 
subspace orthogonal to spanjcfc, • • • , sd}- We denote by {ho, . . . , hk-i) the orthogonal projection 
of (eo, • • • , Gk-i) onto L^. In other words, the hi are defined to be the vectors of M-^+^ which satisfy 
the properties 

hi ■ bj = 6ij Vj < k, and hi ■ Sj = Vj > k (A7) 

One then has: 
Proposition 1 The following equalities hold 

V(eo,--- ,ez)) = V(/io,--- , /ifc-i, e^, • • • ,6^) 

= V(/io, • • • , hk-i) V(efe, • • • , ez)) 

where the function V is defined to be V{vi, ■ ■ ■ , Vn) = det{vi, ■ ■ ■ 

The first equality arises from elementary column operations on the determinant, the second one by 
orthogonality, and the third one from the duality of the basis (/ii)o<i<A;-i and (foi)o<i<fc-i within 
Lk- 

If one applies this result for k = 1 one sees that 

\bi\ = \hi\~^ 
where \hi\ is the height of the D — 1 face i, defined as 

V^ = Z^, (AlO) 

\hi\ 

V being the L)-simplex 'volume' and Vj the volume of the face i. The proposition for k = 2 leads 
to the identity 

VVij = ViVj sin 9ij (All) 

where Vij is the 'volume' of the D — 2 face obtained by dropping i,j and 6ij the dihedral angle 
between the faces i and j. 

Using the well known relation between the cofactor matrix and the inverse matrix one easily 
sees that the dual vectors can be expressed through the derivative of the volume: 

dV 
n = 9^- (A12) 

31 



This implies the derivative formula 



d cos lij 



-2ViVjCos%. (A13) 



This is shown if one looks at the variation 
5V^ = 2V%l5e^i = 2V^ ^(6* • V)ej ■ da = 2V^ ^(6* • V)6{ej ■ a) = -2ViVj ^ cos % 6 cos kj 

i,j i<j i<j 

where summations over repeated indices is understood when it is not explicit, and where i and n 
indices are raised and lowered with flat metric. From (|A1HIAT3]1 . one can compute the derivative 

To show this, lets start from the square of (jAlip and derive it with respect to cos lij. Since neither 
Vi,Vj nor Vij depends on lij, one obtains 

^^' .V2 =V,2v2^i^ = (V,V,sin%)(2V.V,cos%)- ^^'^ 



d cos lij ^^ * ^ d cos lij J J J J d cos lij 

Using once again the identities (|AlllA13p . the RHS is equal to 

dV^ dOii VU, / av^ \ fdOii 



-VV, 



^■' d cos lij d cos lij sin lij y d cos kj J y dkj 



■Kj " ^^3 / ^^ \ / ^^1-0 



from which (JAM]) follows. Finally we derive a last derivative identity by taking the derivative of 
(JAISP with respect to cos lij. With the help of ()A14l lAlip one obtains 



{dcoskj^ ~ ""'^ 



2Vl. (A15) 

^3' 



Next we will show the following 



Proposition 2 We consider a spherical 4-simplex (eo,--- ,64). We denote by lij its lengths and 
G = G(eo, • • • ,64) its Gram determinant. Let Vj be the volume and €j the orientation of the 
tetrahedron Tj obtained by dropping the vertex j. Then: 



dG 



dlQA 
dG 



^^04 



= =F2sin/o4VoV4 (A16) 

'04 

9 V1V2V3 .,,7^ 

= -2eie2e3^ — (A17) 

,± sm to4 

'04 



where {l^^^ei) are solution of G = and uJq^ = mod 27r, with Iq^ < /^. 

The first equality simply arises from (JAISP and from the fact that the 4D dihedral angle of the 
face [123] is vr if /o4 = ^04 and if /o4 = ^04- 

To show the second equality, we first use the Proposition 1 for D = 3 and k = 2 to express the 
volume Vj of the tetrahedron tj, i = 1,2,3, as 

Vi = V(/ij, /ifc)V(eo, 64) = \hj\\hk\ sin 9ij sin Iq4 (A-18) 
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where ijk is a permutation of 123, and \hi\ is the height of the point I in the face [04^]. Indeed, 
first, recall that bj ■ b^ = —\bj\\bk\ cosOjk, where 6jk is the dihedral angle between the faces [04j] 
and [04A:] in r^; second, writing hj ■ h^ = \hj\\hk\ cos 0^ •, the duality of the basis {hj, hk) and {bj, bk) 
within the plane span(6j,6^) implies the relation 0-- = 6ij. 

Now, still with the Proposition 1, the volume of the 4-simplex factorizes as 

V = V(/ii,/i2,/i3)V(eo,e4) = |/ii||/i2||/i3|V(/ii,/i2,/i3) sin/04 (A19) 

Since variations of the dihedral angles leave heights \hi\ and length /o4 invariant, one only has to 
compute the variation of the quantity 

V^(/ii,/i2,/i3) = det(cos6'jj) 

A first way to perform this variation is to use (IA13P with the following correspondence lij — > Oij, 
and Oij -^ aij, where aij is the dihedral angle between edges i and j in the spherical triangle defined 
by {hi,h2,hs)- On can convince oneself that, when G = 0, this spherical triangle is degenerate 
and, furthermore, cos aij = —eiEj. Consequently, one gets 

SV'^iliiMM) = -2eie2e3 sin 012 sin % sin 023 (ei^23 + £2^13 + 63^12) (A20) 

where the term in parenthesis is the variation of the dihedral angle uJq^. 

An other way to obtain this variation is a direct computation of the determinant. Indeed, 

6V (/ii, /i2, /i3) = — 2sin0i2(cos6'i2 — COS013 cos6'23)50i2 + perm (A-21) 

Now, since the deficit angle Uq^ = ei023 + £2^13 + £3^12 vanishes modulo 27r, one can write 

cos 012 — cos 013 cos 023 = COS £3012 — COS 013 COS 023 

= cos(ei023 + £2^13) - COS 013 COS 023 = £162 sin 023 sin 013, (A22) 

and insert this identity in (1A21I) . to get (lA20p . 
Eventually, the variation of G = V^ reads 

5G = —2eie2e2,\hi\ |/i2| I/13I sin /o4 sin 0i2 sin 0i3 sin 023 5a;o4 (A23) 

which, thanks to (jAlSp . can be written as 

Vi V? V3 
6G = -2eie2e3 . ] ^ duj^^ (A24) 

sm to4 

This last equality finally shows ()A17p . 

APPENDIX B: FLATNESS IDENTITIES FOR THE TETRAHEDRON 

We consider a flat tetrahedron IJKL with edge lengths Uj. We define: 0/^ G [0)^] to be the 
geometric angle between the edges {IJ) and {IK) within the triangle IJK, cpjj £ [0, vr] to be the 
exterior dihedral angle of the edge {IJ). We introduce the matrices a/^ = agjK and hfj = /i^e , 

hjj = /i27r-</,=^. , with 

/e*0/2 \ / cos 0/2 sin 0/2 \ 

n6= [ ^ „_,„i/2 ] , ag = I - . - . 1 



e-*<^/V V -sin 0/2 cos 0/2 y 
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Note that, by construction, hij = hjj and a/ = aj . Any SU{2) matrices acts as a rotation on 
M^ by adjoint action g ■ (x'iTj) = g{x'^ai)g~^; as such h^ is a rotation of axis z and angle ^, ag a 
rotation of axis y and angle 9. 

Then for all even permutation (IJKL) of (0123), the following vertex identities hold: 

h^,jaj^'h\Lat''h\^a^J = -el (Bl) 

Moreover, if one considers (jBip as an equation in ^jj where ^jj £ [0,27r] and hjj = h^jj, then 
CiJ = 0/j (the exterior dihedral angle) is the only solution for e = +1, and similarly ^jj = 2Tr — <f>^j 
is the only solution in [0, 2tt] of (jBip for e = — 1. 

To prove this, let us focus on the vertex 0, and let us first give a geometrical interpretation to 
(jBip . We start from a configuration of the edges meeting at such that the edge (01) is along the 
axis z, the edge (02) is in the plane (x, z) and in the direction of the axis x, and the edge (03) is 
in the direction of the axis y. After the rotation a^^ , the edge (02) is along the axis z; after the 
rotation hQ2 the edge (03) is in the plane {x,z). The combination /io2''^o^ J^^t amounts to a cyclic 
permutation of the edges (Oi): after doing three permutations we are back where we have started. 
Therefore the product ()B1|) is equal to ±1. 

In order to establish our statement, lets compute the two sides of the equality: 



In detail the equations are 



ho3afho2 = -e(aJ3)-i(/ioi)-i(42)-i (B2) 



K5oi+5o2+eo3)/2pQg^32/2 = cos0(JV2cos0(JV2-e'««isin0^V2sine^V2 (a) 
<foi+fo2-«o3)/2g-j^032/2 = sin 9^-^/2 cos el'^/2 + 6-'^''^ cos el^/2 sin el'^/2 (6) 



as well as their complex conjugate. The equations (a) and (5) provide four real equations whose 
only three are independent, because of the unitarity condition for the matrices in ()B2p . Indeed, 
considering aa + ab, where a denotes the complex conjugate of the identity (a), leads to the trivial 
equality 1 = 1. Now by taking aa — ab, we get 

C32 = C12C13 - COS^0lSl3'Sl2 (B3) 

and two similar equations obtained by permutation of indices. Cij,Sij denote cos ^q-* , sin ^q-' . These 
equations in S,oi are those satisfied by the exterior dihedral angles. Consequently 

cos% = cos0gj. (B4) 

Therefore, ^oi = (J)q,- or S,oi = 27r — (pQ-, since by hypothesis S,oi € [0, 27r]. 

In order to choose between these two possibilities we need an additional geometrical interpreta- 
tion of these angles: if we surround the vertex by a 2(i-sphere, the three edges meeting at will 
intersect this sphere along three vertex of a spherical triangle, the interior of this triangle being 
inside the tetrahedra. 6q is then the spherical length of the edge (ij) of this triangle, and (/>qj is 
the exterior dihedral angle at the vertex i. The area of this spherical triangle is given by 

A' = 27:-^'oi + <Po2 + ^03, 0<A<2tt. (B5) 

We can take the imaginary part of (a), which implies 

singoi _ cosg|V2 _ 

sin A/2 sin 00^3/2 sin ^0172 " ""^ ^ "^ ^""""^ 
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where ^ = 27r — ^oi + '^02 + Cos- By symmetry on can get similar identities leading to positive 
constants Cj, i = 1 ... 3, and therefore we find 

sinfn,- C,- . , , , 

— |5l = ^>0 V(zj) (B7) 

This last property imposes all the ^oi f o coincide to (/>qj or all to 2tt — (J)q- . In the first case we have 
A = A^\ in the second A = — 27r — A^] in both cases sin A/2 > 0. Therefore one sees that the sign 
of sin^oi is e. This shows that the solutions are (/>oi = (t>Qi when e = +1 and 0oj = 27r — (/)q^ when 
e = -l. 

As an application of the calculation of ()B2p . let us show the expression (j9ip for the jacobian. 
By considering the Euler decomposition of the matrix aQ^/io2fl(}^ = h^agh^ we write 

giCoi) = hoih^aeh^hos (B8) 

The Euler parameters do not depend on ^oi ) Co2 • Therefore 

dg = ^ sin ed0d(eo3 + 0)d(eo2 + ^P) (B9) 

Now we know from previous calculations that the angles ^oi and 9 are related to each other in the 
following way 

cos e = C12C23 - COS Coi'Sl2S23 (Bll) 

where Cij , Sij denotes cos ^q-' , sin 9^^ . Consequently 

sin^d^ = 27rsinCoi'Si2Si3-;^ — (B12) 

= 2.^fl (B13) 

t01'02t03 ^71^ 

V{lij) is (3! times) the volume of the tetrahedron whose edge lengths are kj. The last formula 
holds as long as ^oi is a solution of (jBip . since this angle is then equal, up to a sign and modulo 
2tt, to the exterior dihedral angle of the edge (01). Hence, 

d. = vr-^nT^ (B14) 
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